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Abstract 

We present an optical model (OM) analysis of the elastic scattering data of the re- 
actions ^He+^^Al and ^He+^''*Pb at incident energies around the Coulomb barrier. 
The bare part of the optical potential is constructed microscopically by means of a 
double folding procedure, using the Sao Paulo prescription without any renormal- 
ization. This bare interaction is supplemented with a Coulomb dipole polarization 
(CDF) potential, which takes into account the effect of the dipole Coulomb interac- 
tion. For this CDF potential, we use an analytical formula derived from the semiclas- 
sical theory of Coulomb excitation. The rest of the optical potential is parametrized 
in terms of Woods-Saxon shapes. In the ^He+^^^Fb case, the analysis confirms the 
presence of long range components, in agreement with previous works. Four-body 
Continuum-Discretized Coupled-Channels calculations have been performed in or- 
der to better understand the features of the optical potentials found in the OM 
analysis. This study searches to elucidate some aspects of the optical potential of 
weakly bound systems, such as the dispersion relation and the long range (attractive 
and absorptive) mechanisms. 

Key words: Nuclear reaction ^'^^Fb(^He,^He), halo nucleus. Coulomb dipole 
polarizability. Coulomb barrier, optical potential, threshold anomaly. 
PACS: 25.60.Dz,25.60.Gc,25.60.Bx,21.10.Gv,27.20.+n. 



1 Introduction 



Previous optical model (OM) analyses of elastic scattering data for many 
heavy-ion systems have shown a rapid and localized variation of the optical 
potential at energies in the vicinity of the Coulomb barrier, known as the 
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"threshold anomaly" [T]. This effect is characterized by a rapid increase of 
the surface strength of the imaginary part of the potential with increasing 
energy, up to an approximately constant value, accompanied by a bell-shaped 
peak in the surface strength of the real part. These energy-dependent surface 
strengths are linked by a dispersion relation Physically, this variation 
can be understood in terms of a dynamic polarization potential arising from 
channel coupling effects p!|2|3|^f5] . As the bombarding energy increases, the 
absorption of flux tends to be constant and the energy dependence of the real 
potential is found to be mainly connected with non-local effects [6J. A realistic 
non-local model for the heavy-ion real potential developed in Refs. p[7f5|51[TU] 
has shown to be a powerful tool for studying elastic scattering. In Ref . [10] , a 
simple local approximation to this non-local model was derived. This potential 
factorizes in a fundamental energy- dependent term and an energy-independent 
double folded potential Vp^r), as following: 

VsPp{r) = VF{r)exp[ -], (1) 



where c is the speed of the light and v is the relative velocity between the two 
nuclei. This microscopic model, called Sao Paulo Potential (SPP), has been 
tested for a large number of systems Pl71l5f5puilllfl2fl3fl4] . demonstrating 
its validity at low (sub-Coulomb) and intermediate energies. In particular, it 
has been applied to the scattering of ^He on ^^C and ^^Ni targets, at different 
energies p3|IT6] . without the need of any renormalization of the real bare 
potential. The use of this microscopic potential tends to reduce the ambiguities 
in the real part of the bare potential [T7] . 

Previous works, e.g. Refs. [Ip8] , have shown that the "threshold anomaly" be- 
haviour of the interaction potential involving stable heavy- ion projectiles can 
be explained by coupled-channels calculations employing double folded real 
potentials without renormalizations, provided that the physically significant 
couphngs are explicitly included. These calculations usually employ an inte- 
rior imaginary potential of Woods-Saxon (WS) form to simulate the in-going 
wave boundary condition for fusion, all surface absorption being provided by 
the channels that are explicitly included in the coupling scheme. 

Although the appearance of the threshold anomaly has been studied exten- 
sively in the scattering of light weakly bound ^'^Li [19|20|21ll22|23|24|25j . there 
are very few works done with ^He [26|27] . The ^He nucleus has a weakly bound 
three-body n-n-a Borromean structure and it is known to have an extended 
two neutron distribution. This peculiar structure affects the dynamics of the 
collision, making these processes very interesting, as demonstrated by recent 
experimental campaigns. The data from these experiments show a sizable sen- 
sitivity to this underlying exotic structure [27|28f29|[30] . Reactions induced by 
^He on several targets, at energies around the Coulomb barrier, exhibit some 
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common features, such as a remarkably large cross section for the production 
of a particles. This effect is clearly associated with the weak binding of the 
halo neutrons (S'2n = 0.975 MeV), that favours the dissociation of the ^He 
projectile in the nuclear and Coulomb fields of the target. 

Moreover, we have recently learned that simple preconceptions based on the 
experience of the OM on stable nuclei, such as the role of the strong absorption 
radius, cannot be simply extrapolated to the scattering of exotic nuclei. We 
have also learned about the long range reaction mechanisms which occur when 
^He is scattered on a heavy target, such as ^°^Pb [30|31f32] . These recent 
results confirm that the elastic scattering induced by exotic nuclei on heavy 
targets can be qualitatively different from the scattering of stable nuclei. It is 
observed that the elastic data are sensitive to the values of the real and the 
imaginary potentials in different radial ranges. In general, the sensitivity to 
the real potential corresponds to the region of the strong absorption radius, 
which is Rsa = 12.5 fm for the ^He+^^'^Pb reaction |30]. However, it is found 
that the long range absorption present in the ^He scattering data can only 
be reproduced using large values for the diffuseness of the imaginary part of 
the phenomenological potential, and the sensitivity to the imaginary potential 
occurs at much larger radii [32j. Again, we are able to assess that these data 
can not be well reproduced by short range potentials. 

In the present OM analysis, we try to understand the distinctive features of 
the ^He optical potential, for both light and heavy targets. For this purpose, 
we consider the targets ^^Al and ^°^Pb. In order to keep the underlying ingre- 
dients as fundamental as possible, trying to avoid any kind of ambiguities in 
the real part of the bare interaction, we adopt the double-folding potential of 
Eq. dl]), without any renormalization. This bare interaction is supplemented 
with the analytical Coulomb dipole polarization (CDP) potential of Ref. [53] , 
which describes the effect of the dipole Coulomb interaction on the elastic 
scattering. This CDP potential is entirely determined by the B{E1) distribu- 
tion and does not contain any free adjustable parameter. The remaining part 
of the interaction is described by a complex WS potential, whose parameters 
are determined by a best fit procedure of the elastic data. In order to obtain 
good fits of the ^He-|-^'^^Pb data, both the real and imaginary parts of the 
WS potential require a large value of the diffuseness parameter. In order to 
get further insight of this result, four-body Continuum-Discretized Coupled- 
Channels (CDCC) calculations [34j have been also performed. The trivially 
equivalent local polarization potential derived from these calculations show 
long range attractive and absorptive tails for the real and imaginary compo- 
nents, respectively. This result justifies the need of long-range WS potentials 
in the phenomenological OM analysis of this system. Under this interesting 
scenario, we checked the consistency of the resulting optical potential prescrip- 
tion with the dispersion relation. 
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This paper is structured as follows. In section 2, we first present a brief sum- 
mary of previous ^He,^Li+^''Al and ^He,^Li+^°^Pb OM analyses, pointing out 
some differences between these systems. Afterwards, we present a conventional 
OM analysis based on a complex optical potential, whose real bare potential 
is described by the microscopic SPP potential of Eq. ([1]). This double-folding 
potential is eventually supplemented with an analytical parameter-free, com- 
plex CDF potential which takes into account the effect of dipole Coulomb 
couplings. Unlike the ^B.e+^'^Al case, in order to reproduce satisfactorily the 
®He+^°^Pb data one needs to include real and imaginary long range com- 
ponents in the optical potential. In section 3, we present four-body CDCC 
calculations. The local equivalent polarization potential, derived from these 
calculations, is compared with the phenomenological potentials extracted in 
the OM analysis. Finally, in section 4 we present the summary and conclusions 
of this paper. 



2 Optical Model analysis 

2.1 Previous Optical Model analyses for ^ He reactions 

Within the optical model (OM) framework one must determine appropriate 
values for the geometry of the optical potential. Here we review some previous 
OM analyses of ^He,^Li-|-^''Al,^°'^Pb systems, with the aim of pointing out some 
common and different features of these systems, besides defining our starting 
point. 

It is well known that the ^Li and ^He systems have similar structures. They are 
both weakly bound, with separation energies of 1.475 MeV and 0.975 MeV, 
respectively. The main qualitative difference between them is that the dipole 
Coulomb force can break up ^He into ^He+2n, but it cannot break up ^Li into 
^He-|-^H. The dipole Coulomb operator, in a N = Z nucleus, is an isospin 
1 operator. Since ^Li, ^He and ^H, have isospin in their ground states, it 
is not possible that the dipole Coulomb force breaks up ^Li into ^He+^H. 
Therefore, the mechanisms of ®Li and ^He break up in an intense electric field 
are very different. Breakup of ^Li is governed by the nuclear interaction while 
breakup of ^He is governed by both the nuclear and Coulomb couplings to the 
continuum [5U] . 

Kakuee et al. [22] showed that OM potentials with parameters adjusted to fit 
the ^Li+^°^Pb elastic scattering data are not able to reproduce the ^He-|-^°*Pb 
data at 27 MeV. By contrast, Benjamim et al. [35] have shown that OM cal- 
culations based on the double folding potential of Eq. ([1]) describe satisfactory 
®He-|-^^Al and ^Li+^^Al elastic data at laboratory energies from £^iab=7.0 to 
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a-i (tm) 
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Kei. 


6He+208pb 
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1.75 


1.5 


m 


6He+208pb 


22.0 


5.5(3) 


1.89 


1.9 


m 


6Li+208pb 


109.5 


22.4 


0.88 




[36] 



Table 1 

Real and imaginary WS potential parameters used to fit ^He+^'^^Pb elastic scat- 
tering cross section at 27 MeV from Refs. |31|29j and ^Li+^'^'^Pb elastic scattering 
from Ref. [36] are shown. The rest of parameters are fixed at = 0.811 fm and 
= i^Q. = 7.856 fm. 

13.4 MeV. 

Sanchez-Benitez et al. [31] and Kakuee et al. [32] have performed phenomeno- 
logical OM calculations for ^He on ^ospb at ^iab=14, 16, 18, 22, and 27 MeV, 
using a standard Woods-Saxon (WS) shape. The most remarkable feature of 
the derived potentials is the large value of the imaginary diffuseness parameter 
required to reproduce the data (oj ^ 2 fm). This result was interpreted as an 
evidence of the presence of long-range reaction mechanisms. This is in contrast 
to the ^Li+^°^Pb case, where this phenomenon has not been observed. 

2.2 Optical model analysis for ^He+'^'^^Pb and ^He+'^'^Al 

We start with a conventional OM analysis, using a complex optical poten- 
tial, based on the microscopic Sao Paulo potential (SPP), given by Eq. (1). 
We have considered the experimental data from Refs. [3T|32] for ^He+^°^Pb 
and Ref. [35] for ^He-|-^^A1. The energy-independent part of the bare nuclear 
interaction, Vpir), was calculated with a double folding procedure using the 
effective nucleon-nucleon interaction of Ref. ^lOj, and the matter densities of 
the ^He and ^°*Pb, taken from [T6| and [10], respectively. Both densities are 
parametrized in terms of Fermi-Dirac distributions. For the imaginary part, 
we take the same geometry as the real part, with a normalization factor of 
Nj = 0.78. This is the model OM 1. In the past, this prescription has been 
able to describe a large variety of systems in a very wide range of energies (see 
Ref. [n] for details). 

In Figs. [1] and [2], we compare these OM calculations (solid lines) with the 
experimental data (open circles) for ®He+^^Al and ^He+^°'^Pb systems, re- 
spectively, at several bombarding energies. It can be readily noted that these 
calculations reproduce very well the ^He+^^Al data, in agreement with the 
results reported by Benjamim et al. [35], but they clearly fail to reproduce the 
^He+^°*Pb data. For the higher energies, the calculation shows a rainbow peak 
that is not present in the data. Moreover, the elastic scattering cross section 
is overestimated by the calculations, mainly at backward angles, suggesting 
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Fig. 1. Elastic scattering angular distribution in the center of mass frame for the 
reaction ^He+^^Al at different bombarding energies (£'iab=9.5, 11, 12, and 13.4 
MeV). The circles are the experimental data from Ref. [35]. The OM calculation, 
using double-folding real and imaginary components (OM 1), are displayed by solid 
lines. The dot-dashed lines are the four-body CDCC calculations. By dashed lines 
are shown the OM calculations using the OM 1 model plus the CDP potential. 
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Table 2 

Results of the OM analysis for the reaction ^He+^'^^Pb using an energy-independent 
geometry. For all the fits, the only free parameter is the depth of the imaginary 
potential, wl- The radius and the diffuseness of the imaginary potential were fixed 
to Rio = 7.856 fm and Oj = 1.9 fm, respectively. There is no renormalization in the 
bare potential. 

that the flux removed from the elastic scattering due to the reaction channels 
is underestimated. 

To improve the agreement with the data, we have performed a second analy- 
sis, in which we keep the double-folding SPP for the real part of the nuclear 
interaction, but we describe the imaginary part by means of a standard Woods- 
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Fig. 2. Elastic scattering angular distribution in the center of mass frame for the 
reaction ^He+'^'^^Pb at different bombarding energies (i^iab=14, 16, 18, 22, and 27 
MeV). The circles are the experimental data from Refs. |32|31j . The different OM 
model analyses are presented by solid (OM 1), dashed (OM 2) and dotted (OM 3) 
lines. The four-body CDCC calculations are shown by dot-dashed lines. 
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Table 3 

Results of the OM analysis for 6He-H^°^Pb at energies from 14 to 27 MeV, using 
an energy-dependent geometry for the phenomenological imaginary potential. The 
only free parameters are the depth and diffuseness of the imaginary potential. The 
radius of the imaginary potential was kept fixed to RiQ = 7.856 fm. There is no 
renormalization in the bare potential. 

Saxon (WS) shape, with adjustable parameters. We also include an interior 
imaginary potential with a WS shape to simulate the in-going boundary con- 
dition for fusion. So, in this new analysis the projectile-target potential is 
parametrized as: 

U{r) = VsMr) + ^Ws{r) + ziy^r), (2) 
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with 



Ws{r) 



1 + exp 



1 + exp (^) 



(3) 



The parameters of the interior WS potential {Ws{t)) are kept fixed at all 
energies to the values: ws = 50 MeV, Rg = 3.87 fm, and = 0.2 fm. 

The parameters of the external WS potential {WL{r)) are adjusted for each 
energy to fit the elastic data. With the aim of checking consistency, we as- 
sumed the same geometry of the imaginary WS potential obtained by Sanchez- 
Bem'tez et al. [3l] (Table ^ for the ^He+^ospb system, obtained at 22 and 27 
MeV, where we expect the elastic scattering to be most sensitive to the ge- 
ometry of the nuclear potential. Thus, the reduced radius and the diffuseness 
of the imaginary potential were fixed to Rio = 7.856 fm and a, = 1.9 fm, 
respectively. 

With the optical potential geometry fixed, we allowed the imaginary depth 
(wl) to vary in order to reproduce the elastic scattering data at laboratory 
energies of 14, 16, 18, 22, and 27 MeV. These OM fits have been performed 
with the routine SFRESCO, which is part of the FRESCO code §7\. This 
is the model OM 2. The extracted values of the depths are summarized in 
Table [2] and the corresponding angular distributions are compared with the 
data in Fig. |2] (dashed lines). The fits reproduce very well the experimental 
angular distributions, but with some limitations at 27 MeV and backward 
angles at 22 MeV. 

As a third step, and in order to study the sensitivity of the data to the ge- 
ometry, we have allowed the diffuseness of the imaginary potential to vary, 
along with the depth. The diffuseness value remained around 1.9 fm (1.6 - 
2.6 fm). The extracted parameters are summarized in Table |3l This is model 
OM 3. The new fits, shown in Fig. |2]by dotted lines, provide a better agree- 
ment with data. However, at 22 MeV, this calculation shows a pronounced 
rainbow, which does not follow the trend of the data with the energy at the 
measured angles. These large values of the imaginary diffuseness parameter 
are consistent with the prescription of Bonaccorso and Carstoiu [38] . Using a 
semiclassical approach, and assuming that the imaginary part of the optical 
potential is well represented by an exponential form, they estimated that the 
diffuseness of the imaginary potential is approximately given by Oj ~ (27j)~^, 



with 7j = y2/i|eb|, where Eb = Ti^^f is the binding energy. If this formula 

is applied to ^He, one gets ^ 2 fm, which is in very good agreement with 
the results of the present work. 

The importance of this diffuse imaginary potential is better seen in Fig. |3l 
where we plot the strength of the imaginary WS component as a function of the 
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Fig. 3. Strengths of different imaginary optical potentials that fit the data, as a 
function of the interacting distance, for the ^He+^'^^Pb system. The diffuseness 
parameter for the imaginary part used in the calculations are in the vicinity of the 
minimum obtained from the best fit analysis (a, = 1.50 — 3.0 fm). 

interacting distance, for different values of the diffuseness parameter a^. The 
values of at chosen for this plot are close to the minimum (a^ = 1.5 — 3 fm). 
Despite the well known ambiguity in the choice of the diffuseness parameter 
[T7] we find that, in order to reproduce satisfactorily the data set, one has to 
use a large value of the diffuseness parameter a,. In addition, for each energy, 
the radius of sensitivity of this potential, defined at the distance at which the 
different potentials cross, is well beyond the strong absorption radius, meaning 
that the elastic scattering will be affected by the details of this potential at 
large distances. 

2.3 Optical Model analysis including a Coulomb dipole polarization potential 

It is well known that the elastic scattering of weakly bound nuclei on heavy 
targets is strongly affected by the polarization induced by the dipole part 
of the Coulomb interaction. In particular, it has been found that this dipole 
polarizability effect gives rise to a significant reduction of the elastic scattering 
cross sections, which is particularly important for the collision of weakly bound 
nuclei on heavy targets. This effect has been shown to account for part of the 
long range behaviour found in phenomenological OM analyses of these systems 
[31.32]. 
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Here, we investigate how this effect varies for two extreme cases of the charge 
of the target, and in which case we identify long range Coulomb polarization 
as a distinctive feature of the scattering of ^He, at energies around the bar- 
rier. Within the OM framework, the effect of the dipole polarizability on the 
elastic observables can be included by means of a Coulomb dipole polarization 
(CDP) potential. A simple analytical expression for this CDP potential was 
derived in [33|39] . The form of the polarization potential is obtained in a semi- 
classical framework requiring that the second order amplitude for the dipole 
excitation-deexcitation process and the first order amplitude associated with 
the polarization potential are equal for all classical trajectories corresponding 
to a given scattering energy. This leads to an analytic formula for the polar- 
ization potential for a single excited state [33]. The expression so obtained can 
be generalized for the case of excitation to a continuum of breakup states [39] 
giving rise to the following formula: 



2^2 



Upoi{r) 



9 hv {r — a^Yr 



(4) 



X de 



dB{El,e) 
de 



£6 



where oq is the distance of closest approach in a head-on collision, v is the 
projectile velocity and g and / are analytic functions defined as 



f{z, = ie^' exp i-nOK^,^ ^z) , (5) 

— oo 

and ^ = ^ is the Coulomb adiabaticity parameter corresponding to the ex- 
citation energy e of the nucleus. An important feature of this potential is 
that when the breakup energy Sb is large enough, the purely real adiabatic 
dipole potential is re-obtained. In the opposite limit, for small breakup ener- 
gies f (- — l,n — > 1 and g (- — — )■ 0, and the polarization potential 
becomes purely imaginary, depending on r as ^a^ - For our analysis we 
take theoretical values of the B{E1) distribution of ®He [40ll4T] . This deter- 
mines completely the CDP potential at different energies without using free 
parameters. 

In this section, we reanalyze the ^He+^^Al and ^He-l-^'^^Pb elastic data, includ- 
ing explicitly the effect of dipole polarizability by means of the CDP potential 
described above. In the ^He+^^Al case, the optical potential contains the dou- 
ble folding SPP potential, an imaginary part with the same geometry, and the 
CDP potential. The real and imaginary parts of the double-folding potential 
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Fig. 4. Elastic scattering angular distribution in the center of mass frame for the 
reaction ^He+'^'^^Pb at different bombarding energies (£'iab=14, 16, 18, 22, and 27 
MeV). The circles are the experimental data from Ref. |32|31j . The different OM 
model analyses, using the prescription of Eq. ([7]), are presented by solid (OM 4) 
and dashed (OM 5) lines. 
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Table 4 

Best-fit optical potential parameters for the system ^He+^''^Pb, including explic- 
itly the CDF potential and using energy-independent geometries for the optical 
potential. For all the fits, the only free parameter is the depth of the imaginary po- 
tential. The radius and diffuseness parameter of the imaginary potential were fixed 
to RiQ = 7.856 fm and aj = 1.1 fm, respectively. 

were renormalized by Nr = 1.0 and Nj = 0.78, respectively. The results of 
these calculations are represented in Fig. [1] by dashed lines. As expected, we 
find a negligible effect of the CDP potential in the elastic scattering angular 
distributions. 



For the ^He+^°^Pb reaction, the optical potential contains the real double- 
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folding SPP potential (with Nr = 1), the CDP potential, and WS imaginary 
components, i.e.: 



(7) 



The long range component Wiir) will account for the effect of other periph- 
eral reaction channels not included in the CDP potential. As before, we have 
divided our procedure into two steps. In the first search, we allowed the imagi- 
nary geometry to vary for each energy. In a second search, we try to keep fixed 
the optical potential geometry with the energy, searching for the minimum. 
The results of these calculations are shown in Fig. |H The solid lines (OM 4) 
correspond to the calculations keeping fixed the geometry of the WS potential, 
whereas the dashed lines (OM 5) are the calculations allowing the diffuseness 
parameter to vary. In both cases, we observe a similar agreement with data. 
The depths for the optical potentials extracted with the energy-independent 
geometry are listed in Table HI Comparing these results with those of the pre- 
vious subsection, we see that the inclusion of the CDP potential produces a 
similar quality of data fits, but leads to a reduction of the imaginary diffuse- 
ness from 1.9 fm to about 1.1 fm. We can conclude that the long range tail 
found in the OM analysis of the ^He-|-^°'^Pb data is partially due to the effect 
of dipole Coulomb breakup. Although the overall agreement with the data is 
good, it seems that some ingredient is missing, since we still have some limi- 
tations to reproduce the data at energies above the Coulomb barrier (Fig. H] 
at 22 and 27 MeV). 

2.4 Optical Model analysis with a long range real potential 

Given the impossibility to reproduce the data above the barrier, we propose in 
this subsection a more general prescription in which, besides the phenomeno- 
logical imaginary part, we introduce also a phenomenological real part of WS 
shape. In these new serial of calculations, we keep the complex analytical CDP 
potential to include Coulomb breakup effects. Consequently, this analysis was 
performed using the following parametrization for the optical potential: 



where Ws{r) and Wiir) describe, as before, the short range and long range 
components of the imaginary nuclear polarization potential, and VL{r) is a 
real phenomenological polarization potential. 



U{r) = Vspp{r) + ?7cDp(r) + iWs{r) + VL{r) + ^W^i(r), 
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Fig. 5. Elastic scattering angular distribution in the center of mass frame for the 
reaction ^He+'^'^^Pb at different bombarding energies (£'iab=14, 16, 18, 22, and 27 
MeV). The circles are the experimental data from Ref. |32|31j . The different OM 
model analyses are presented by solid (OM 6) and dashed (OM 7) lines. 

As in our previous calculations, we do not renormalize the bare potential. The 
strengths vl and w^, as well as the diffuseness parameters of the WS potentials 
(or and Oj) have been allowed to vary (OM 7). Afterwards, we searched the 
values of and Oj that best fit the angular distributions at all energies, and we 
fixed them (OM 6). The best data fits are obtained using the values = 1.1 fm 
and Oj = 1.0 fm. The extracted parameters are summarized in Tableland the 
corresponding angular distributions are presented in Fig. OM 6 by solid lines 
and OM 7 by dashed lines. It can be seen that these calculations reproduce 
very well the data in the whole angular and energy range, with no significant 
differences between the two calculations. This analysis suggests that, besides 
the inelastic couplings produced by the dipole Coulomb force, there are other 
long range mechanisms that give rise to a diffuse tail in the phenomenological 
optical potential. 

These calculations show a complex scenario from where we could expect impor- 
tant changes in the dispersion relation of this system as compared to normal 
systems. In Fig. El we plot the strength of the energy-independent geometry 
optical potential (as a total sum of the real and imaginary OM components), 
as a function of the collision energy for the ^He+^'^^Pb system, and for different 
values of the projectile-target distance. The results seem to show certain cor- 
relation in the variation of the real and imaginary parts, almost independent 
on the interacting distance, which is qualitatively consistent with dispersion 
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relations [2]. However, the errors bars limit our conclusions. These errors bars 
correspond to a variation on the complex long range potential around the Xmin 
that results on an increase of the total by an amount Xmin/-^- This analysis 
claim for better precision experiments in order to infer about the existence of 
a threshold anomaly in reactions involving the weakly bound ^He nucleus. 

The variation of the real and imaginary potentials with the energy has been 
studied also in terms of their volume integrals. For this purpose, in Table [5] 
we include also the values of the real and imaginary volume integrals ( and 
J^) per interacting nucleon pair, defined as: 



where Ap, At are the projectile and target masses. and Uyj are, respec- 
tively, the total real part (excluding the monopole Coulomb contribution) and 
the total imaginary part (excluding the interior Woods-Saxon potential) of 
the optical potential. For the real part, the volume integral remains roughly 
constant, slightly decreasing with increasing incident energy. This energy de- 
pendence is consistent with existing parametrizations (see eg. |l2|43|44j ) and 
is a consequence of Passatore's |15|16] application of Feshbach dispersion re- 
lation. Within the energy interval considered in the present work, the energy 
variation of is small, remaining around the value 360 MeV fm^, which is 
consistent with the value found by Mohr [17] in the analysis of the ^He+^°^Bi 
reaction at energies around the Coulomb barrier. By contrast, the values of 
the imaginary volume integrals show a more pronounced variation with the 
bombarding energy. According to the dispersion relations |48j, the decrease of 
Jy should be accompanied by an increase of J^. In the present case, there is 
not a clear trend in the energy dependence of the extracted values of J^. Nev- 
ertheless, one has to keep in mind that this energy dependence comes mainly 
from the Woods-Saxon potentials Vl and Wl and, as we have shown above, 
the strengths extracted for these potentials have large error bars. 

The large errors bars shown in Fig. [6] for the real and imaginary strengths of 
the phenomenological potentials might lead to the conclusion that the features 
found by these potentials could be just a consequence of the ambiguities of 
the real and imaginary parts. However, in the next section, we show that 
these features, namely, the long range behaviour of the real and imaginary 
components, arise from physical effects due to specific reaction mechanisms 
present in the collision. 
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Table 5 



Best-fit parameters for the reaction °He+ Pb, including explicitly the CDP po- 
tential and the real and imaginary WS potentials VL{r) and WL^r). For all the fits, 
we considered an energy-independent geometry, with Rq = RiQ = 7.856 fm, a,- = 1.1 
fm, and Oj = 1.0 fm, which provide the best overall fit. For each energy, the free 
parameters are the depths vl and wl- The columns Jy and refer to the real and 
imaginary volume integrals, as described in the text. 
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Fig. 6. Real and imaginary strengths (V, W) of the optical potential as a function 
of the bombarding energy, for the system ^He+^'^'^Pb. The strengths are evaluated 
at the distance indicated by the labels. The real potential strength is the sum 
of the double folding potential Eq. ([T]), the real part of the CDP potential and 
the phenomenological component Viir). The strength of the imaginary potential is 
the sum of the WS potentials of Eq. ([3]) and the imaginary part of the dynamic 
polarization potential. 
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3 Microscopic analysis within a few-body approach 

3.1 Continuum-Discretized Coupled- Channels calculations 

Theoretically, the treatment of the reactions induced by the Borromean nu- 
cleus ^He requires a four-body formalism (three-body projectile plus a target). 
For three-body problems (two-body projectile plus a target) the Continuum- 
Discretized Coupled- Channels (CDCC) framework [l^fSU] has been success- 
fully used for many cases [5T|52] . Recently this method has been extended 
to four-body problems [53|54f55|l56] . In particular it has been tested for the 
reaction ^He on ^"^^Pb at 22 MeV [31]. Consequently, in this work we use this 
four-body CDCC framework for performing theoretical calculations for the 
reactions under study, ^He+^°^Pb and ^He+^^Al at different energies around 
the Coulomb barrier. We have used the binning procedure [31] for discretizing 
the three-body continuum of the projectile. This procedure has been devel- 
oped very recently as an extension of the same method used traditionally in 
standard three-body CDCC calculations. For this purpose, the three-body con- 
tinuum representation uses the eigenchannel expansion of the multi-channel 
S-matrix [3^ . 

Here we use the same structure model for the three-body system ^He(a-|-n+n), 
as in Refs. [5HI5B] . The Hamiltonian includes two-body potentials plus an ef- 
fective three-body potential. Continuum states with angular momentum and 
parity j'^ = 0"*",!", and were considered. The wavefunctions for these states 
were generated using the codes FaCE [57J and sturmxx [58|. The maximum 
hypermomentum used was -ft'max = 8. The parameters of the three-body in- 
teraction are adjusted to reproduce the ground-state separation energy and 
matter radius (for j = 0^ states) and the resonance energy (for j = 1~ and 2+ 
states). The calculated ground state energy was 0.953 MeV and the root mean 
squared (rms) radius was 2.46 fm (assuming a rms radius of 1.47 fm for the a 
particle). Both Coulomb and nuclear potentials are included. The fragment- 
target interactions were represented by optical potentials which reproduce the 
elastic scattering at the appropriate energy. The n+^^^Ph and n+^^Al poten- 
tials were from [52], the a-l-^^^Pb was from from [20], and for a+^^Al we used 
the code by S. Kailas [61], which provides optical model parameters for a 
particles using real potential volume integrals of Atzrott et al. ^8], geometry 
systematics, and dispersion relation. 

The coupled-channels equations were solved using the code FRESCO [57] . 
that reads the coupling potentials externally. We included in the calculation 
the projectile-target interaction multipole couplings with order Q = 0,1, 2. In 
order to get convergence, the number of eigenchannels included was 4 for both 
systems. However, the maximum energy value £max, the number of bins nbin for 



16 



System 


^lab (MeV) 


£max (MeV) 


nbin (0+,l-,2+) 




Rrn (fm) 


6He+208pb 


14.0 


5.0 


(6,9,6) 


150 


200 


6He+208pb 


16.0 


6.0 


(9,12,9) 


150 


200 


6He+208pb 


18.0 


7.0 


(9,12,9) 


150 


200 


6He+208pb 


22.0 


8.0 


(6,9,6) 


150 


200 


6He+208pb 


27.0 


8.0 


(6,9,6) 


150 


200 


6He+27Al 


9.5 


5.0 


(6,9,6) 


30 


80 


6He+27Al 


11.0 


6.0 


(6,9,6) 


30 


80 


6He+27Al 


12.0 


7.0 


(6,9,6) 


30 


80 


6He+27Al 


13.4 


8.0 


(6,9,6) 


30 


80 



Table 6 



Parameters for the four-body CDCC calculations. See text for details. 

each j'^, the maximum total angular momentum Jmax and the matching radius 
Rm depended on the target and on the energy. These values are presented in 
Table d 

Figs. [1] and [2] show the four-body CDCC calculations for the ^He+^^Al and 
^He-f ^°*Pb systems, at different energies around the Coulomb barrier (dash- 
dotted lines). For the ^He+^''A1 case, the calculations reproduce well the trend 
of the data, although some underestimation is observed. For the ^He+^'^^Pb 
case, these calculations reproduce very well the forward angular region. How- 
ever, as the energy decreases the calculation tends to underestimate the data 
at larger angles. For this reaction, we have used a second set of optical poten- 
tials for the fragment-target interactions. These were generated by the SPP 
using Eq. ([T]). The calculations with these potentials yield very similar results, 
although there are small differences at backward angles. However the magni- 
tude of these differences cannot explain the underestimation of the data at 
these large angles. 

3.2 Local equivalent polarization potential from CDCC calculations 

In order to link the results of this section with those found in the OM analysis, 
we have extracted from the four-body CDCC calculations the so called triv- 
ially equivalent local polarization (TELP) potential [IB]. This is a local and L- 
independent potential which represents the overall effect of the breakup chan- 
nels on the elastic scattering. This potential is constructed in such a way that 
the one-channel calculation performed with the potential f/bare('^) + ^telp('^) 
gives the same elastic scattering as the full CDCC calculation. The bare po- 
tential, t/bare('^) IS just the sum of the fragment-target interactions convoluted 
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Fig. 7. ^He+2°^Pb (left) and ^He+^'^'Al (right) TELP potentials extracted from the 
four-body CDCC calculations. 

with the ground state density of the ^He nucleus. Figure [7] shows these polar- 
ization potentials (due to Coulomb and nuclear interaction) calculated for the 
different systems at different energies. 



For the ^He+^^Al system (right panel in Fig. [7]), we observe a very uniform 
behavior of the dynamic polarization potential (real and imaginary parts) as a 
function of the interacting distance. The real part of the dynamic polarization 
potential is repulsive (except at very short distances where the details of this 
potential are probably not meaningful). Moreover, this potential shows a weak 
dependence with the incident energy. This result explains why the use of a 
double folding bare potential and an imaginary part with the same geometry 
is able to reproduce satisfactorily the data. 

On the contrary, for the ^He+^°^Pb system (left panel in Fig. [7]), such uni- 
form behavior is not observed. The real part is repulsive at distances close to 
the strong absorption radius and becomes attractive at larger distances. The 
imaginary part is mostly absorptive, although for the lowest energies it be- 
comes emissive at short distances!!] Both the real and imaginary parts extend 
to large distances, well beyond the strong absorption radius. These features 



It has been pointed out [62] that this emissive imaginary part is a consequence of 
representing a strongly non-local object, namely, the dynamic polarization potential 
arising from the coupled-channels couplings, by a simple local potential. This effect, 
nevertheless, does not lead to unitary breaking. 
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are consistent with the findings of Mackintosh and Keeley [62] and Rusek [63] 
for the same reaction. Clearly, this complicated behaviour cannot be simply 
simulated by a renormalization of the double-folding potential. 

To get a deeper understanding on the relationship between the features of the 
TELP and the phenomenological potentials extracted in the OM analysis of 
the ^He+^°^Pb reaction we show in Fig. |8] the separate contributions of the 
TELP potentials arising from either Coulomb (left) or nuclear (right) cou- 
plings. It can be seen that the Coulomb couplings are responsible for the long 
range attractive tail in the real part of the TELP potential. Also, these cou- 
plings produce a long range absorptive tail. Within the phenomenological OM 
analysis, this behaviour is expected to be at least partially taken into account 
by means of the analytical CDP potential. Also from Fig. [8] one sees that nu- 
clear couplings are responsible for the strong repulsive part of the polarization 
potential. Although they are of shorter range than the Coulomb polarization 
potential, it is noticeable that both the real and imaginary components of 
the nuclear polarization potential extend also to distances well beyond the 
strong absorption radius. Therefore, besides long range Coulomb couplings, 
the ^He+^°^Pb reaction is characterized by long range nuclear couplings. In 
particular, the long range absorptive tail explains why the OM calculations 
using the double-folding potential for the imaginary part do not reproduce 
the data (even after inclusion of the CDP potential), and supports the need 
of a long range complex component in the OM potential. This long range 
part of the real polarization potential cannot be well accounted for by a mere 
renormalization of the double-folding potential. 

It is worth to note that the conclusions extracted from the TELP have to 
be analysed with caution, since this potential is just a local L-independent 
approximation of a very complicated non-local and L-dependent object. In 
particular, the oscillations of the TELP at short distances arise from the ra- 
dial solution of the elastic equation, which appears in the denominator in the 
expression of the TELP. To support the conclusions extracted from this poten- 
tial, we have extracted another polarization potential fitting the elastic angular 
distribution obtained with the CDCC method with a local optical potential. 
This potential contains the double-folding SPP potential and Woods-Saxon 
real and imaginary components, with fixed radius R=7.86 fm. The parame- 
ters obtained from this fit {vl=-32.1 MeV, 0^=0.98 fm, wl=10.2 MeV, ai=1.7 
fm) confirm that, in order to reproduce the effect of the continuum couplings 
using a single-channel optical potential, the real and imaginary parts need a 
long-range absorptive tail. 

It is illustrative to compare these results with those of previous works. The 
effect of the breakup channels on the optical potential has been subject of 
many studies in the past. Although it is not the aim of this paper to make an 
extensive review of these works, we cite some previous results which are closely 
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Fig. 8. ^He+^°®Pb Coulomb (left) and nuclear (right) TELP potentials extracted 
from the four-body CDCC calculation. 



related to ours. In the eighties, the Kyushu group studied in detail the effect of 
the continuum in the elastic and transfer of deuterons using the CDCC method 
[iniEDj- By extracting the polarization potential from the CDCC calculations, 
they concluded that the breakup induces a surface complex polarization po- 
tential with a real repulsive part and an imaginary absorptive part. This is 
consistent with our results for the light target. However, for the lead target, 
we have a significant long-range attractive component and the imaginary part 
is of much longer range than that found in the cited work. As we have already 
pointed out, this is a consequence of the strong Coulomb couplings. These 
couplings where omitted in these pioneering calculations by the Kyushu work 
and, in addition, the effect is expected to be much smaller for the deuteron 
case due to its larger binding energy. 

Similar conclusions where achieved in the comprehensive work of Sakuragi 
[6l] for ^Li scattering. He found that the breakup channels produce a strong 
repulsive term and an absorptive part. The values of these potentials in the 
strong absorption radius were consistent with the renormalization required 
to reproduce the elastic data using double-folding potentials. His conclusions 
can not be readily extrapolated to the present case because (i) only nuclear 
breakup was included and (ii) in any case the effect of dipole couplings are 
very much suppressed in ^Li with respect to ^He. Nevertheless, his results are 
consistent with our calculations when Coulomb breakup is switched off. 

Finally, we cite the work of Matsumoto et al. [M] where they study the elastic 
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scattering of ^He on ^°^Bi at Coulomb barrier energies within a pseudo-state 
version of the four-body CDCC method. They find the real part of the equiv- 
alent local polarization potential, extracted from their CDCC calculations, 
is repulsive at short distances and becomes attractive for distances beyond 
^15 fm, whereas the imaginary part is absorptive and of long range. These 
results are in total agreement with ours, as expected given the similitude be- 
tween both reactions. 



4 Summary and conclusions 

Detailed optical model (OM) analyses of the ^He on ^^Al and ^ospj,, ^jg^i^g^^ 
at laboratory energies around the Coulomb barrier (i?iab=9.5, 11.0, 12.0, and 
13.4 MeV and ^iab=14, 16, 18, 22, and 27 MeV, respectively), have been 
performed. 

In the case of the light system, ^He-|-^^A1, we cannot recognize breakup ef- 
fects on the optical potential analysis. Thus, a conventional optical potential, 
with the same form factor for the real and imaginary parts, based on the 
double-folding potential given by Eq. ([T]), have been successfully applied to 
describe the data of this light weakly bound system. The inclusion of a com- 
plex Coulomb dipole polarization (CDP) potential, which takes into account 
the effect of the dipole Coulomb interaction, had no significant effect on the 
data fits. Four-body Continuum-Discretized Coupled- Channels (CDCC) cal- 
culations corroborate the conclusions of the OM analysis. 

In the case of the heavy system, ^He-|-^°^Pb, breakup effects (Coulomb and 
nuclear) are very relevant and play an important role in the dynamics of 
the reaction around the Coulomb barrier. These couplings give rise to long 
range attractive and absorptive components in the optical potential required 
to reproduce the data. Despite this complex scenario, we tried to develop an 
optical potential, with an energy-independent geometry, consistent with the 
dispersion relation and suitable for studying reaction mechanisms. 

For the analysis of the ^He-|-^°^Pb data, we adopt the microscopic real nuclear 
Sao Paulo potential given by Eq. ([1]), with its energy- independent geometry 
calculated as a double-folding of the effective nucleon-nucleon interaction [10] 
with the ®He and ^osp]^ matter densities. We performed calculations with and 
without including explicitly the analytical CDP potential. Both the double- 
folding and the CDP components have no adjustable parameters. The fact 
that, even after the inclusion of the CDP potential, one still requires relatively 
intense values of the optical potential at large distances, means that there are 
other relevant long range mechanisms, besides the dipole Coulomb polariz- 
ability. To obtain a satisfactory agreement with the data above the barrier, a 
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repulsive long range Woods-Saxon (WS) potential (with real and imaginary 
parts) had also to be included in the OM analysis. In this approach, the op- 
timal values for the diffuseness parameters were = 1.1 fm and = 1.0 fm. 
These values are to be compared, respectively, with the diffuseness value of 
the nuclear densities of the double-folding potential, which is = 0.56 fm, 
using a Fermi-Dirac representation, and with the mean imaginary diffuseness 
obtained to fit the data set, without using the CDP potential, Oj = 1.9 fm. 
In particular, this reduction of the imaginary diffuseness when including the 
CDP potential confirms that Coulomb breakup accounts for a considerable 
part of the long range behavior found in the OM analysis. 

The existence of sizable long range effects in the ^He+^'^^Pb interaction, sug- 
gests the presence of reaction mechanisms that remove flux from the elastic 
channel at distances well beyond the strong absorption radius {Rsa = 12.5 fm) 
and even well below the Coulomb barrier. These results can easily be under- 
stood since, as it has been discussed in this work and in previous ones (see, 
e.g. Ref. [30]), in the case of ®He+^°^Pb reaction, the breakup cross section 
has its maximum at 18 MeV and decreases for energies around (below and 
above) this value, but it remains very large for these energies. It evidences 
the presence of a strong dynamic polarization potential, which has shown to 
be consistent with the trivially equivalent local polarization (TELP) potential 
derived from four-body CDCC calculations (Figs. [7] and [8]). For projectile- 
target separations around the strong absorption radius, the real part of the 
TELP potential is very repulsive (Fig. [7]), while it becomes attractive at large 
distances (Figs. [7]and[H]). The repulsive part has been identified as coming 
from nuclear couplings, whereas the attractive part arises from Coulomb cou- 
plings. The presence of the repulsive part in the polarization potential justifies 
the inclusion of the real WS potential in the optical model analysis. Both nu- 
clear and Coulomb couplings produce also a long range absorptive component 
in the TELP potential, which again justifies the inclusion of the imaginary 
WS component in the optical model calculations. Moreover, the complicate 
behaviour of the TELP potential as a function of the distance, as well as its 
strong energy dependence, suggest that this complicated behaviour cannot be 
simply accounted for by a renormalization of the double folding potential. 

Therefore, our analysis suggests a new and complex scenario for the optical 
potential involving the scattering of the weakly bound ^He nucleus on heavy 
targets. The attractive and absorptive effects, produced by the couplings to the 
continuum breakup states, tend to produce significant changes in the strength 
of the real and imaginary parts of the optical potential with the energy, at the 
different interacting distances. Consequently, any OM analysis involving this 
nucleus as well as the application of the dispersion relation to the study of 
threshold anomaly, must be made with caution. With the aim of assessing the 
consistency with the dispersion relation, depending mainly on the dynamic 
polarization potential effects, we have studied the energy dependence of the 
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real and imaginary parts of the potential evaluated at different values of the 
interacting distance. The results, presented in Fig. |6l show a certain correla- 
tion in the variation of the real and imaginary parts with the energy, almost 
independent on the interacting distance, which is consistent with the disper- 
sion relation [2]. However, our conclusions are limited by the errors bars and 
so, in order to allow us to infer about the dispersion relation and threshold 
anomaly in the ^He reactions, more accurate data would be required. 

For the future, couplings between the transfer/breakup channels and the elas- 
tic channel must be incorporated beyond the first order, thus performing a 
coupled-reaction channel calculation. This calculation could allow an assess- 
ment on whether the explicit inclusion of these channels can account for the 
remaining part of the long range absorption effects. 
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